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Two energy scales of spin dimers in clinoclase Cu3(As04)(OH)3 
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Magnetic susceptibility and microscopic magnetic model of the mineral clinoclase 
Cu3(As04)(OH)3 are reported. This material can be well described as a combination of 
two nonequivalent spin dimers with the sizable magnetic couplings of J ^ 700 K and Jd2 — 300 K. 
Based on density functional theory calculations, we pinpoint the location of dimers in the crystal 
structure. Surprisingly, the largest coupling operates between the structural CU2O6 dimers. We 
investigate magnetostructural correlations in CU2O6 structural dimers, by considering the influence 
of the hydrogen position on the magnetic coupling. Additionally, we establish the hydrogen 
positions that were not known so far and analyze the pattern of hydrogen bonding. 
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I. INTRODUCTION 

The majority of magnetic systems have more than one 
characteristic energy scale, according to the different na- 
ture of interactions between the spins. For example, 
in low-dimensional magnets strong interactions within 
chains or planes are of direct exchange or superexchange 
type, whereas weak interchain (interplane) couplings may 
have purely dipolar origin.^ ^ The different energy scale 
of these interactions implies that, at sufficiently high tem- 
peratures, the magnetic behavior is solely determined by 
the strong couplings, and the system can be fully de- 
scribed in terms of a low-dimensional spin model. ^ ^ How- 
ever, at low temperatures, interchain (interplane) cou- 
plings come into play, and a full three-dimensional de- 
scription is required. 

In complex spin systems, the identification of different 
energy scales is by no means a simple problem. Naively, 
one could think that distinct crystallographic positions of 
magnetic atoms should lead to different strengths of mag- 
netic couplings and, therefore, to different energy scales 
of the magnetic behavior. Indeed, in spin-dimer systems, 
such as BaCuSi206 and NH4CUCI3, the spin-^ Cu^+ 
ions occupy several nonequivalent positions^ ^ and form 
different types of spin dimers which have large impact 
on the high-field behavior, including the unusual criti- 
cal regime of the Bose-Einstein condensation of magnons 
in BaCuSi2 06 (Refs.|7j) and the fractional magnetization 
plateaus in NH4CUCI3 (Refs. "S" and "9). However, spin 
dimers can also be formed between two nonequivalent 
Cu positions, thus leading to only one type of spin dimer 
and one energy scale, as in the spin-ladder compound 
BiCusPOeP 

When the system contains several Cu^+ positions with 
dissimilar local environment and variable connectivity of 
the Cu polyhedra, the identification of relevant interac- 
tions and energy scales becomes increasingly complex. 
The problem of magnetic dimers that do not match struc- 
tural dimers, ^^^^ as well as magnetic chains running 
perpendicular to the structural chains ,1^ is well-known 
in quantum magnets and requires a careful microscopic 



analysis. 

Here, we report on the magnetic behavior and micro- 
scopic modeling of clinoclase, Cu3(As 04)(O H)3. The in- 
tricate crystal structure of this mineraP^^ features three 
nonequivalent Cu positions. Nevertheless, the resulting 
spin lattice comprises only two types of magnetic dimers 
with notably different interaction energies. Our micro- 
scopic analysis shows that one spin dimer is formed be- 
tween two different Cu positions (Cul-Cu2) and does 
not match the respective structural dimer. However, 
the other spin dimer coincides with the Cu3-Cu3 struc- 
tural dimer. We argue that neither a straightforward 
comparison of Cu-Cu distances nor the application of 
the Goodenough-Kanamori-Anderson (GKA) rules for 
the superexchange^^ ^^ lead to the correct assignment 
of magnetic couplings in clinoclase. Therefore, all geo- 
metrical details of relevant exchange pathways should be 
taken into account and considered simultaneously. We 
elaborate on this problem by determining the positions 
of hydrogen atoms and analyzing their role in the su- 
perexchange. 

The outline of the paper is as follows. Sec. |TI| sum- 
marizes experimental and computational techniques that 
were applied in this study. In Sec. |III A] we report details 
of the sample characterization followed by the determi- 
nation of hydrogen positions and the discussion of the 



crystal structure in Sec. IIIB Sec. IIIC presents experi- 
mental magnetic properties and their brief discussion in 
terms of a phenomenological model of two spin dimers 
with different energy scales, whereas Sec. |III D] provides 
a microscopic insight into this model and into residual 



interactions between the spin dimers. Finally, Sec. HIE 
clarifies the role of hydrogen atoms in the Cu-O-Cu su- 
perexchange. Our work is concluded with a discussion 



and summary in Sec. IV 



II. METHODS 

For our experimental studies, we used a natural sample 
of clinoclase from Wheal Gorland, St. Day United Mines 
(Cornwall, UK), which is the type locality of this rare 



mineral. The sample provided by the mineralogical col- 
lection of Salzburg University (Department of materials 
engineering and physics) features bulky dark-blue crys- 
tals of clinoclase together with smaller light-blue crys- 
tals of liroconite, Cu2Al(As04)(OH)4-4(H20). The crys- 
tals of clinoclase were manually separated from foreign 
phases and carefully analyzed by powder x-ray diffrac- 
tion (XRD) and chemical analysis. Laboratory powder 
x-ray diffraction (XRD) data were collected using the Hu- 
ber G670 Guinier camera (CuKq, i radiation, ImagePlate 
detector, 20 = 3 — 100° angle range). Additionally, high- 
resolution XRD data were collected at the ID31 beamline 
of the European Synchrotron Radiation Facility (ESRF, 
Grenoble) at a wavelength of about 0.43 A. The chemical 
composition was determined by the ICP-OES method. ^^ 
The thermal stability of clinoclase was investigated by 
thermogravimetric analysis^ ^ up to 500 °C. 

Magnetic susceptibility of clinoclase was measured 
with a Quantum Design MPMS SQUID magnetometer 
in the temperature range of 2 — 380 K in fields up to 5 T. 

Electronic structure calculations were performed 
within density functional theory (DFT) by using the full- 
potential local-orbital code FPL09. 07-41.^^ Local density 
approximation (LDA)^and generalized gradient approx- 
imation (GGAp^were used for the exchange-correlation 
potential, together with a well converged /c-mesh of 
6x6x6 points for the crystallographic unit cell of clin- 
oclase and about 100 points in supercells. Hydrogen po- 
sitions missing in the presently available crystallographic 
data^^ ^^ were obtained by structural optimizations with 
a threshold for residual forces of 0.002 eV/A. 

The effects of strong electronic correlations, typical for 
cuprates, were considered by mapping the LDA bands 
onto an effective tight-binding (TB) model. The trans- 
fer integrals ti of the TB-model are evaluated as non- 
diagonal matrix elements in the basis of Wannier func- 
tions (WFs). These transfer integrals ti are further in- 
troduced into the half-filled one-orbital Hubbard model 
H = Htb + K^ff X^i ^it^ii ^hat is eventually reduced to 
the Heisenberg model for the low-energy excitations. 



(ij) 



(1) 



where the summation is done over bonds {ij). For the 
half-filled case, which applies to clinoclase, the reduction 
to the Heisenberg model is well-justified in the strongly 
correlated limit ti <C t/eff, with the effective on-site 
Coulomb repulsion [/eff exceeding ti by at least an or- 
der of magnitude (see Table |ll| . This procedure yields 
AFM contributions to the exchange coupling evaluated 
as J^FM = ^tyjj^^^ 

Alternatively, full exchange couplings J^, comprising 
FM and AFM contributions, can be derived from differ- 
ences in total energies of various collinear spin arrange- 
ments, as evaluated in spin-polarized supercell calcula- 
tions within density functional theory LSDA+/7 formal- 
ism that includes a mean-field correction for correlation 
effects. An "around mean field" (AMF) as well as a 



"fully localized limit" (FLL) approximation for correct- 
ing the double counting were used.^^ Both types sup- 
plied consistent results so that only the AMF results are 
presented here. The on-site Coulomb repulsion and on- 
site Hund's exchange of the Cu Zd orbitals are chosen as 
Ud = 6.5±0.5 eV and Jd = ^ eV, respectively, according to 
the parameter set used for several other cuprates. ^^^^^ 

In addition to periodic DFT calculations, we per- 
formed a series of cluster calculations that pinpoint the 
effect of hydrogen atoms on the superexchange. The clus- 
ter under consideration is based on the CU2O6H5 dimer 
and embedded in a set of point charges, with two As 
ions bonded t o the dimer considered as total ion po- 
tentials (TIPs) pHEol r^Yie embedding was chosen so that 
the intradimer hopping obtained from the cluster and 
periodic LDA calculations match. The cluster calcula- 
tions were done with the Orca 2.9 cod d^^^^^ in com- 
bination with a 6-311-^+G(d,p) basis set and a PBEO 
hybrid functional.'^ 

Quantum Monte Carlo (QMC) and simulations were 
performed using the codes LOOP^ and DlRLOOP_SSEp^ 
from the software package ALPS-LS.!^ Magnetic sus- 
ceptibility was simulated on finite lattices of A^=2400 
spins 5* = I , using 40 000 sweeps for thermalization and 
400 000 sweeps after thermalization. For simulations of 
magnetization isotherms, we used 4000 sweeps for ther- 
malization and 40 000 sweeps after thermalization. Mag- 
netization of the "2+1" model was simulated using the 
full diagonalization code from ALPS-l.S.'^ 



III. RESULTS 
A. Sample characterization 

Powder XRD confirmed the purity of our clinoclase 
sample. However, the ICP-OES analysis showed slight 
deviations from the ideal composition: 49.6(3) wt.% 
Cu, 18.3(1) wt.% As compared to 50.1(1) wt.% Cu, 
19.7(1) wt.% As expected for Cu3(As04)(OH)3. Addi- 
tionally, trace amounts of Ca and S (0.1 — 0.2 wt.%) were 
found. Other detectable elements, including transition 
metals, are below 0.03 wt.%. 

To verify the composition of the clinoclase phase, we 
performed structure refinement using high-quality syn- 
chrotron data collected at room temperature. The lattice 
parameters we obtained for the space group P2i/c are 
a = 7.266 A, 6 = 6.459 A, c = 12.393 A with the mono- 
clinic angle /3 = 99.49°. They agree well with the existing 
single-crystal data,^^ where a = 7.257 A, 6 = 6.457 A, 
c = 12.378 A and (3 = 99.51°, which we used for the mag- 
netic modeling. Our refinement confirmed full occupa- 
tion of all atomic positions in the clinoclase structure.^ 
As Cu and As are the two heaviest elements (and, there- 
fore, strongest scatterers) in Cu3(As04)(OH)3, their con- 
tent in the crystalline phase is safely established by XRD. 
Regarding the bulk composition, the slight deficiency of 
Cu and As, as revealed by chemical analysis, may be 




FIG. 1. (Color online) The crystal structure of clinoclase. The Dl (Cul-Cu2) structural dimers are shown in orange, the D2 
(Cu3-Cu3) dimers in greenish color. The left panel shows a single layer consisting of zig-zag chains of Dl-dimers linked by 
As04 tetrahedra (pink). In the right panel, the sandwich- like structure of clinoclase is visible with the hydrogen bonds shown 
as (blue) dashed lines. Arrows indicate the leading hopping pathways. The central panel shows high-quality natural clinoclase 
crystals from the Majuba Hill Mine, Pershing Co., Nevada, USA. 



attributed to trace amounts of secondary phases, such 
as CaCOa and CaS04 that are possible impurities in 
natural samples. Another plausible impurity is CuO 
(tenorite), which is difficult to identify by XRD, because 
its strongest reflections overlap with those of clinoclase. 
Note that none of the possible impurities reveal any con- 
spicuous effects in the magnetic susceptibility and should 
not affect our experimental results reported in Sec. [Ill C) 
Thermogravimetric analysis (TGA) identified the on- 
set of the decomposition at about 180 °C and the weight 
loss of 7.1 % upon heating to 500 °CW^ This weight loss 
implies the release of 1.5 water molecules per formula 
unit, as expected for Cu3(As04)(OH)3. The sample re- 
covered after the heating contained a mixture of unknown 
phases. Their composition and crystal structures require 
further investigation that lies beyond the scope of the 
present study. 



B. Crystal structure 

Clinoclase crystallizes in the monoclinic space group 
P2i/c and forms a fairly complex crystal structure with 
three nonequivalent Cu positions (Fig. [I]).^^ All Cu sites 
have a five- fold square-pyramidal coordination, where the 
four oxygen atoms in the basal plane form shorter bonds 
to Cu (1.9 — 2.1 A), whereas the distance to the axial oxy- 
gen atom is above 2.3 A. Similar to other cuprates, su- 
perexchange pathways and ensuing magnetic interactions 
can be described by resorting to the planar CUO4 coor- 
dination, because the oxygen atom in the axial position 
of the pyramid does not take part in the superexchange 
(see Sec. [mD ). 

The CUO4 plaquettes around the Cul and Cu2 sites 
form doubly bridged dimers (Dl = Cul-Cu2) that share 
corners and build zigzag chains directed along [001]. The 
ASO4 tetrahedra connect these chains into layers in the 
be plane. Two layers of this type form a "sandwich" 



encompassing the dimers of Cu3 atoms (D2 = Cu3-Cu3). 
Although both Dl and D2 dimers are built of two CUO4 
plaquettes sharing a common edge, their symmetries and 
geometrical parameters are different. For example, the 
inversion symmetry in the center of D2 entails two equal 
Cu3-0-Cu3 angles. In contrast, Dl has no symmetry 
elements, hence its two Cul-0-Cu2 bridging angles are 
not equal. 

The connections between the "sandwiches" are re- 
stricted to hydrogen bonds and to long Cu-0 bonds in 
the axial positions of the CuOs pyramids. This weak in- 
terlayer bonding is responsible for the perfect cleavage 
of clinoclase crystals parallel to (100), and for the low 
(Mohs)-hardness of 2.5-3 in this material.'^ 

Previous structure determinations^^ were based on 
XRD data and did not report the positions of hydrogen 
atoms. However, precise positions of all atoms, including 
hydrogen, are required for DFT band structure calcula- 
tions and ensuing microscopic analysis of the electronic 
and magnetic structure. Therefore, we determined the 
hydrogen positions by relaxing the crystal structure of 
clinoclase. Only the hydrogen positions were optimized, 
whereas all other atoms were fixed to their experimen- 
tal positions.l^S In the starting model, hydrogens were 
attached to three out of seven oxygen atoms at a typi- 
cal 0-H distance of 1.0 A. The orientation of the 0-H 
bonds was random, although we made sure that the hy- 
drogens are well separated from other atoms in the clin- 
oclase structure. While there is freedom in choosing three 
oxygen atoms forming covalent bonds to hydrogen, only 
those oxygens that do not belong to the ASO4 tetrahe- 
dra led to structures with low energies. When hydrogens 
are attached to oxygens belonging to the ASO4 tetrahe- 
dra, the energy is much higher, hence such structures 
can be ruled out. This is in agreement with the empir- 
ical assignment of the OH groups in the experimental 
crystallographic study.l^ 

The resulting hydrogen positions are listed in Table [l[ 



TABLE I. Hydrogen positions obtained by LDA/GGA struc- 
ture optimization. The positions of oxygen forming short 0-H 
bonds to these hydrogen atoms are given in brackets. 



atom 



x/a 



y/h 



z/c 



HI (05) 0.7429/0.7493 
H2 (06) 0.9362/0.9346 
H3 (07) 0.1522/0.1511 



0.3513/0.3515 
0.4685/0.4680 
0.1635/0.1648 



0.4799/0.4807 
0.6792/0.6786 
0.4876/0.4864 



Further details of the relaxation procedure and compar- 
isons to the experiment for other Cu^+ hydroxy-salts 
can be found in Refs. [26] and [40] The optimized 0-H 
distances are close to 1.0 A, as expected for the cova- 
lent 0-H bonds. Each hydrogen atom also forms one 
longer contact of about 1.8 A (hydrogen bond) to an- 
other oxygen atom. Two of these contacts provide addi- 
tional bonding within the layer, ^^ whereas the hydrogen 
bond formed by H2 connects adjacent "sandwiches" (see 
Fig. [l]). This arrangement of hydrogen bonds correlates 
with the positions of the hydrogens: while HI and H3 
are nearly coplanar with Cu and O atoms, H2 is no- 
tably shifted along the a direction toward the neighbor- 
ing "sandwich" . Therefore, HI and H3 lie in the plane of 
the Dl dimer plaquettes, whereas H2 and the respective 
0-H bond are in the out-of-plane position. 



C. Magnetization measurements and 
phenomenological fits 

The temperature dependence of the magnetic suscep- 
tibility x{T) is shown in Fig. [5] In quantum magnets, 
x(T) typically has an asymmetric dome-like shape, with 
a broad maximum indicating a gradual crossover from the 
high-temperature (paramagnetic) to the low-temperature 
(correlated) regime. Due to the unusually high magnetic 
energy scale in clinoclase, this maximum is shifted to high 
temperatures (~300K) and is barely visible in the data 
collected below and around room temperature. Unfortu- 
nately, high-temperature measurements are not possible 
because the decomposition of clinoclase starts at about 
450 K.^^ 

The upturn in x(^) below 50 K is a typical extrinsic 
feature caused by defects and/or impurities. It can be 
reasonably described by a Curie- Weiss law with Cimp = 
0.015 emu K /mol, corresponding to 3.2% oi S = \ im- 
purities per f. u., and 6>imp = 2.5K. After subtraction of 
the extrinsic contribution, we obtain vanishingly small 
susceptibility below 30 K and an activated (exponential) 
behavior at higher temperatures, evidencing the gapped 
nature of the magnetic excitation spectrum. 

The gap between the lowest-lying 5 = and S = 1 
states (the spin gap) is inherent to numerous magnetic 
models. The simplest one is a quantum-mechanical spin 
dimer with a singlet ground state. Indeed, the structure 
of clinoclase features well-defined structural dimers Dl 
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FIG. 2. Experimental magnetic susceptiblity (circles) of 
clinoclase Cu3(As04)(OH)3 and fits using a dimer model 
(dash-dotted line), the phenomenological model [Eq.ls] of 2-hl 
dimers (solid line) , and the microscopic model of coupled spin 
dimers, as shown in Fig. Is] (dashed line, QMC fit). Besides 
the intrinsic dimer susceptibility [Eq. [2] , we accounted for the 
temperature-independent contribution and a Curie- Weiss im- 
purity and/or defect contribution. Inset: difference curves for 
the 2+1 dimers solution and the QMC fit. 



and D2 that are evocative of the spin-dimer magnetism. 
For a system of isolated dimers, the magnetic suscepti- 
bility is given by the exact analytical expression: 



X{T) 



1 



T (3 + exp[J/T])' 



(2) 



where J is the magnetic exchange within the dimer. The 
fit yields J = 415 K, but it does not account for the shape 
of the experimental curve, as shown in Fig. [2] Moreover, 
the resulting value of the ^-factor (^ = 1.48) is unrealisti- 
cally small for Cu^+. 

The failure of the isolated dimer model suggests that 
the magnetic couplings in clinoclase are more intricate. 
Prior to the microscopic evaluation (Sec. HID), we at- 
tempt to describe the x(^) dependence with a phe- 



nomenological magnetic model. According to Sec. IIIB[ 
the monoclinic unit cell of clinoclase contains three struc- 
tural dimers formed by edge-sharing CUO4 plaquettes 
(Fig. IT]). Two of these dimers are of the Dl type (Cul- 
Cu2), whereas the third dimer is D2 (Cu3-Cu3). Natu- 
rally, the magnetic couplings within Dl are different from 
those in D2. This brings us to a tentative model of 2+1 
dimers: 



H, 



Dl 



H = 2Hdi + Hd2 
- 22 ^Di(Si -S^) 



(3a) 
(3b) 



iGCul 

jeCu2 



H, 



D2 



^ Jd2(S.-S,), (3c) 

{iJ)eCu3 

where Jdi and Jd2 denote the magnetic couplings in the 
respective dimers. For each dimer, the magnetic suscep- 
tibility is given by Eq. Q. 



Two scenarios are possible: either Jdi > Jd2, or the 
other way around. Fitting to the experimental curve 
readily shows that the Jd2 > ^di solutions do not repro- 
duce the experimental behavior. In contrast, the model 
with Jdi > Jd2 yields an excellent fit with Jdi = 703.5 K 
Jd2 = 289.3 K, and g = 1.86 (fuh line in Fig. |2|. There- 
fore, the smaller gap of ~290 K comes from the Cu3-Cu3 
dimers (D2), while the Cul-Cu2 dimers (Dl) give rise to 
the larger gap of ~700 K. 



D. Microscopic magnetic model 

Now, we compare the above phenomenological model 
with the microscopic results based on DFT. In a first step, 
LDA calculations are performed. The width of the whole 
valence band block of about 9 eV is typical for cuprates 
(see Fig.[3|. The spurious metallicity of the energy spec- 
trum is a well-known shortcoming of LDA due to the un- 
derestimated electronic correlation in the Cu 3d shell.'^ 
Nevertheless, the LDA bands around the Fermi level (the 
energy range from —0.5 to 0.7 eV) are sufficient to de- 
scribe low-energy magnetic excitations, provided that a 
suitable correlation part is added to the model Hamilto- 
nian. The relevant bands are essentially of the Cu dxy 
character, with sizable contributions from O 2p orbitals. 
The orbital symmetry is defined with respect to the local 
coordinate system on each Cu04-plaquette, where the 
Cu-Cu bond of the dimer is chosen as the local x-axis, 
and the z-axis is orthogonal to the plaquette plane. Note 
that this setting is different from the standard one, where 
X and y axes are directed along the Cu-0 bonds, so that 
the highest crystal- field level has the dx'2-y2 symmetry. 

The leading hopping parameters ti and correspond- 
ing AFM exchanges J/^^^ = 4tf//7eff are listed in Ta- 
ble |TI| The results of the model analysis are supported 
by the evaluation of full exchange integrals Ji using total 
energies of collinear spin configurations calculated with 
LSDA+/7. These two approaches are complementary. 
The model analysis provides information on all exchange 
couplings in the system and guides the LSDA+[/ calcu- 
lations that are restricted to only a handful of leading 
interactions. 

Our model analysis based on the hopping parameters ti 
identifies five leading AFM exchanges that exceed 100 K 
(see J/^^^ in Table M. The perfect agreement between 
the LDA bands at the Fermi level and those calculated 
with the Cu-centered Wannier functions (Fig.jsl confirms 
that the relevant superexchange pathways in clinoclase 
can be well described in terms of the CUO4 plaquettes. 
Despite the short Cu-Cu distance of only 3.3 A between 
the Cu atoms in two contiguous "sandwiches" along the a 
direction, the respective hopping is negligibly small (be- 
low 20 meV) because the magnetic orbitals lie in the be 
plane and do not overlap. Likewise, two outer layers of 
each "sandwich" are coupled only via the Cu3 spins and 
lack any direct interaction. 

The comparison between jf^^^ and Ji in Table O 
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FIG. 3. (Color online) The LDA density of states (DOS) 
and the band structure of clinoclase Cu3(As04)(OH)3. The 
top panel shows the contributions of the Cu(3d), 0(2p) and 
As states to the total DOS. The Fermi level is at zero energy. 
In the bottom left panel, the LDA-bands around the Fermi 
level are displayed and compared with bands derived from 
a fit using an effective one-band tight-binding model based 
on Cu-centered Wannier functions (WFs) projected on local 
Cu(3(iccy)-orbitals. The /c-points are defined as follows: F = 
(000), X= (fOO), Y= (OfO), Z= (OOf ), XZ= (JO^),MZ= 
(?f f )' YZ= (Of f ). The bottom right panel shows that the 
partial Cu(3(i)-D0S at the Fermi level is basically of Cu(3dxy) 
character, justifying our construction of the WFs. 



TABLE II. Leading exchange couplings in clinoclase 
Cu3(As04) (011)3: Cu-Cu distances dcu-Cu (in A), bridging 
angles (/?cu-o-Cu (in deg), hopping integrals ti (in meV), an- 
tiferromagnetic contributions j/^^^ = 4t^/[/eff (in K) with 
[/eff = 4.0eV. The Ji (in K) are calulated with LSDA+[/ us- 
ing Ud = 6.5 ± 0.5eV, Jd = leV. 





C?Cu-Cu 


V^Cu-O-Cu 


t^ 


jAFM 


J 


Jdi 


2.98 


93.6/99.9 


-115 


153 


-4 ±8 


Jd2 


3.13 


101.9 


191 


423 


302 ± 53 


Jidi 


3.38 


124.6 


-117 


159 


161 ±25 


J 


3.66 


149.3 


276 


884 


693 ± 99 


Jid2 


5.52 


- 


-106 


130 


159 ±31 



shows that four out of five leading interactions are indeed 
AFM, with only small FM contributions. However, the 
coupling within Dl is nearly canceled because of compa- 
rable FM and AFM terms. The large FM contribution to 
^Di (J™ = Jdi - Jm^ ^ -160 K) is indeed expected 
for the coupling geometry with bridging angles close to 
90°. 

Our microscopic model is consistent with the phe- 
nomenological analysis that suggested spin dimers with 
two different energy scales (Sec. IIIC). The coupling on 



D2 is Jd2 = 302 K very close to 300 K found experi- 
mentally. However, the magnetic dimer with the cou- 
pling of about 700 K has to be re-assigned. The cou- 
pling within Dl is in fact very weak, so that the spin 
dimer with J :^ 700 K is formed not on the Cul-Cu2 
structural dimer Dl, but between the respective dimers, 
where the CUO4 plaquettes share a common corner in- 
stead of sharing a common edge. This effect can be well 
understood in terms of the GKA rules for the superex- 
change, because the Cu-O-Cu angle for J is nearly 150° 
compared to only 93 — 100° for Jdi- However, the GKA 
rules do not account for the fact that Jd2 — 302 K ex- 
ceeds Jidi — 159 K, even though the bridging angle for 
Jd2 is notably smaller. This unusual behavior is further 
discussed in the next sections. 

The re-assignment of the magnetic dimer has no ef- 
fect on the fit of the magnetic susceptibility presented 
in Sec. [HI C] since it is independent of the position of the 
dimers in the crystal structure. The interdimer couplings 
Jidi and Jici2 are non-frustrated and can be taken into ac- 
count by QMC. The resulting fit shown in Fig. [2] is only 
slightly better than the fit with the phenomenological 
"2+1" model. We find J = 706.8 K, Jd2 = 318.1 K and 
g = 1.893 in good agreement with our previous results. 
The coupling between the magnetic dimers were chosen 
as Jidi = Jid2 = 0.125 J to yield best agreement with the 
experimental curve. 



E. Role of hydrogen 

One important difference between the coupling path- 
ways for Jdi, Jd2, and Jidi pertains to the positions of 
hydrogen atoms. While the bridging oxygen atoms of 
D2 belong to the ASO4 tetrahedra and have only weak 
contacts to hydrogen, the bridging atom for Jidi, as well 
as one of the bridging atoms of Dl, form covalent 0-H 
bonds. For studying the role of the 0-H bonds in more 
detail, we will focus on the Dl dimer where, according to 
the very small exchange coupling and the planar CU2O6 
geometry, we expect interesting effects, as, e.g., a change 
from FM to AFM coupling. 

The effect of the out-of-plane angle r of the 0-H bond 
on the intradimer coupling was studied by Ruiz et al)^ 
for small organic ligands. They found that a large r (out- 
of-plane position of hydrogen) favors FM coupling. We 
attempted to verify this effect for the Dl-dimer in clin- 
oclase. In a first step, the ^di hopping parameters are 




FIG. 4. (Color online) Wannier functions (WFs) on the Cul 
(yellow-green) and Cu2 sites (red-blue). The net overlap at 
the two bridging oxygen within the Dl dimers is significantly 
smaller than the overlap at the oxygen bridging two dimers. 
The Cul-WF exhibits a considerable distortion towards the 
As04-tetrahedra responsible for the large Jid2 coupling. 




FIG. 5. (Color online) Wannier functions (WFs) on the Cu2 
(red-blue) and Cu3 (grey- violet). The overlap between Cu2 
and Cu3 is hampered by the non-planar arrangement of the 
Cu04 plaquettes. The Cu3-WF features distortions of the 
0(2p) contributions caused by the ASO4 tetrahedra. 



calculated as a function of r in the periodic model. The 
out-of-plane rotation of H up to 67° , which corresponds 
to the optimized crystal structure of clinoclase, reduces 
^Di by about 40% and, thus, the AFM contribution to 
Jdi by about 60% (see Supplemental material^). Fur- 
thermore, t is slightly reduced and the inter-sandwich 
hoppings decrease by about 50%. All other hoppings are 
more or less independent of r. 

For the calculation of Jdi as a function of r, we used 
the CU2O6H5 cluster model, embedded in TIPs and point 
charges, which allows to investigate the intradimer cou- 
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FIG. 6. (Color online) The effect of the out-of-plane rotation 
of hydrogen on the Wannier functions (WFs) localized at the 
Cul and Cu2 sites. The top figure shows the Dl structural 
dimer with H, bonded to the bridging oxygen, being rotated 
out of the dimer plane by an angle r. The central and bottom 
figures show the WFs at the two Cu sites for r = 7° and r = 
67°, respectively. The Cu04 plaquettes are indicated in light 
blue, black lines connect the central Cu and the O bonded 
to H. With increasing r, the 0{2p) orbital rotates (indicated 
by orange arrows) about a local z-axis, which is visible with 
respect to the black line. Black arrows point at the H(ls) 
contribution to the WFs. At r = 67° that corresponds to the 
computationally relaxed structure, no such contributions are 
visible anymore. 



pling exclusively. Additionally, the cluster enables us to 
vary the bridging angles without changing the whole set 
of additional structural parameters, as this would be the 
case for a periodic model. ^^ 

The results of the cluster calculations (Fig. u\ nicely 
show the transition from AFM to FM coupling upon an 
increase in r. This effect is driven by the reduced hop- 
ping, because the FM contribution J™ = Jbi — ^of ^ 
is weakly dependent on r and hovers around —150 K. 
The absolute size of Jdi obtained in the cluster calcula- 
tions with the PBEO functional is somewhat larger than 
the LSDA+/7 estimates. This is in fact little surpris- 
ing, because hybrid functionals, such as PBEO, tend to 
overestimate the exchange couplings.'^ 

The decrease in toi can be traced back to the in- 
creasing contribution of the bridging oxygen 0{2py) to 
the WFs of the Cul and Cu2 sites (Fig. [t]), while the 
0{2px) contribution remains constant. At small r, H(l5) 
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FIG. 7. (Color online) The upper panel shows Jdi as a 
function of the out-of-plane angle of hydrogen, r. In the lower 
panel, the ratio 0{2py)/0(2px) of the WF-contributions of 
the bridging oxygen bonded to H are displayed. The orbital 
character is denoted with respect to the coordinate system 
shown in Fig. [6] H(t)/H(t = 7) shows the H(ls) contribution 
normalized to the value at r = 7°. The vertical dashed line 
indicates the computationally relaxed out-of-plane angle of 
hydrogen. 
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FIG. 8. Microscopic magnetic model of clinoclase 

Cu3(As04)(OH)3. Note that each structural "sandwich" de- 
picted in Fig. [l] comprises three interpenetrating lattices of 
this type that are decoupled from each other. For the nota- 
tion of magnetic couplings, see Table III] 



strongly interacts with 0{2py) and thus shifts its orbital 
energy downwards, which in turn reduces the interac- 
tion between this oxygen orbital and Cu{3dxy) orbitals. 
The H(ls) contribution itself decreases with increasing 
r. This is also visible in the WF-picture (Fig. p| as the 
rotation of the contribution of bridging oxygen atoms. 
As the r increases, the 0{2p) orbital turns into the di- 
rection perpendicular to the Cu-Cu axis of Dl, hence the 
overlap of the WF's of Cul and Cu2 is reduced. 
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FIG. 9. Magnetization isotherm of clinoclase 

Cu3(As04)(OH)3 simulated using exact diagonalization of the 
"2+1" dimers model and QMC for the microscopic magnetic 
model with J : Jd2 : Jidi : Jid2 = 1 : 0.45 : 0.125 : 0.125. Mag- 
netic field is scaled adopting J and g values from the x(^) 
fits (Fig.[2|. msat is the saturation magnetization. Note two 
wide plateaus at ?7i/??isat=0 and ?7isat = l/3. 



IV. DISCUSSION 

The magnetism of clinoclase is well described by a 
model of two nonequivalent spin dimers. While this 
model is easi ly der ived from the crystal structure of the 
mineral (Sec. IIIC), the arrangement of magnetic dimers 
cannot be established on purely empirical grounds. Two 
shortest Cu-Cu distances are formed within the struc- 
tural dimers Dl and D2. Assuming that the magnetic 
coupling is more efficient at short distances, one would 
identify these structural dimers as magnetic dimers. This 
assessment is correct for D2, yet it fails for Dl, where the 
strong interaction J forming the magnetic dimer is found 
between the structural Dl dimers. The GKA rules could 
provide a more plausible scenario, because they account 
for the fact that short Cu-Cu distances entail Cu-O-Cu 
angles close to 90° that are unfavorable for an AFM cou- 
pling. The reference to the GKA rules readily explains 
why the dimer is formed by J (bridging angle of 149.3°) 
and not by Jdi (bridging angles below 100°). However, 
a consistent application of these rules would also relegate 
Jd2 (101.9°) that should be well below Jidi with a much 
larger Cu-O-Cu angle of 124.6°. 



In clinoclase, neither Cu-Cu distances nor Cu-O-Cu 
angles fully elucidate the microscopic scenario. To ex- 
plain why Jd2 exceeds Jidi^ details of their superex- 
change pathways should be analyzed and compared. 
Apart from the hydrogen atoms considered in Sec. |III E[ 
we find two main differences between these couplings: i) 
the number of bridging oxygen atoms, which is two for 
Jd2 and one for Jidi ; ii) the mutual arrangement of the 
Cu04 plaquettes that are coplanar for Jd2 and strongly 
twisted for Jidi . Dividing the overall Cu-Cu hopping ^02 
by two, we obtain t^JI = 95 meV that reflects the transfer 
via a single Cu-O-Cu bridge (as in Jidi)- This hopping 
is slightly below \tidi\ = 117 meV, but their difference 
is much smaller than expected for the bridging angles of 
101.9° and 124.6°, respectively. For example, the model 
calculation from Ref. ^441 suggests that tidi should be at 



least twice larger than ^02- 

There are different scenarios explaining the large AFM 
coupling of Jd2- It may originate from the combined ef- 
fect of the indirect Cu-O-Cu and direct Cu-Cu hoppings 
within the CU2O6 structural dimer. While the Cu-O-Cu 
processes should be solely determined by the bridging an- 
gle, the direct hopping requires the coplanar arrangement 
of the CUO4 plaquettes, which is the case for Jd2 only. 
This explanation is in line with the robust AFM coupling 
observed in many other spin-dimer compounds, such as 
TICUCI3 (Ref. 45) and SrCu2(B03)2 (Ref. 46), despite 
their low bridging angles of 96 — 98°. Same arguments 
could be applied to Dl, where one of the bridging angles 
is as large as 100° and indeed leads to a sizable trans- 
fer ^Di = —1 15 me V. However, the out-of-plane 0-H 
bond (see Sec. HIE) has strong impact on the superex- 
change and is responsible for the very weak coupling. 
This effect of side groups may also play an important 
role for D2, yet in a different manner. Here, the bridging 
oxygen atoms belong to the ASO4 tetrahedra that could 
amplify the AFM superexchange, similar to Ge04 tetra- 
hedra in CuGeOs.^^ Note that the structural dimers in 
a-Cu2As207 (bridging angle of 101.7°) are ve ry sim ilar 
to D2 and also feature a strong AFM coupling .^^^1^ 

Another interesting feature of clinoclase is the sizable 
interdimer coupling Jid2 • In contrast to all magnetic cou- 
plings discussed so far, Jid2 does not involve a direct con- 
nection between the CUO4 plaquettes and occurs via the 
bridging ASO4 tetrahedron. The efficiency of this su- 
perexchange pathway can be explained by the coplanar 
and, moreover, well aligned arrangement of the CUO4 
plaquettes. Their positions are such that two Cu-0- • • O- 
Cu contacts are formed. One of these contacts goes along 
the edge of the ASO4 tetrahedron, while the second con- 
tact does not involve any bonds or polyhedra. Neverthe- 
less, its short O- • • O distance of 2.8 — 3.0 A is likewise 
beneficial for the superexchange. The resulting coupling 
Jid2 — 140 K is comparable to the typical interaction 
via double PO4 and ASO4 bridges in Sr2Cu(P04)2!^ 
K2CuP207,^and Cu2As207>^ 

Regarding the spin model of clinoclase, the "strong" 
(J) and "weak" (Jd2) magnetic dimers are joined into 
a planar structure (Fig. [8| by non-frustrated interdimer 
couplings Jidi and Jid2' Three interpenetrating planes 
of this type together form one structural "sandwich" and 
remain nearly decoupled. Topologically, each of the three 
lattices represents a diluted square lattice of magnetic 
dimers, depicted schematically in Fig. |8] 

This spin lattice only marginally differs from a sim- 
ple superposition of two nonequivalent dimers. The phe- 
nomenological "2+1" dimer model and the 2D spin lat- 
tice provide nearly indistinguishable fits of the magnetic 
susceptibility (Fig. [2]). Field dependence of the magneti- 
zation, as calculated by QMC, is very close to the intu- 
itive picture of isolated spin dimers (Fig. [9|. The wide 
plateau at M = | is due to the saturation of the "weak" 
dimers {J-d2)^ while the "strong" (J) dimers remain in the 
singlet state. Therefore, Jidi and Jid2 have little effect 



on the magnetic susceptibility of clinoclase. 

The effect of these weak couphngs is visible by com- 
paring the magnetization curves simulated for isolated 
dimers and for the dimers coupled by Jidi and Jid2^ 
as shown in Fig. |8] In the model augmented by Jidi 
and Jid2^ the transitions preceding and following the |- 
plateau are broadened compared to the "2+1" dimer 
model (Fig. [9]). This broadening is caused by the in- 
terdimer couplings Jidi and Jid2 that give rise to a dis- 
persion of magnetic excitations, and mediate magnon- 
magnon interactions underlying the peculiar effect of 
spontaneous magnon decay.^ 

Clinoclase belongs to the family of gapped quantum 
magnets with nonequivalent spin dimers. In contrast 
to other systems of this type, different spin dimers are 
inherent to the crystal structure of this mineral. They 
are not formed due to a symmetry reduction upon a 
low-temperature phase transition that keeps similar mag- 
netic interactions and, therefore, similar spin gaps in all 
dimers, as in BaCuSi206 (Ref.[5]) and NH4CUCI3P Clin- 
oclase can be rather compared to the ambient-pressure 
modification of (VO)2P2 07, where two distinct spin gaps 
of 32 K and 65 K define two different energy scales of 
the system. ^^^^ Systems of this type may show inter- 
esting high-field behavior, because each group of dimers 
(J and Jd2) has independent low-temperature transi- 
tions related to the Bose-Einstein condensation (EEC) 
of magnons. The BEC transition takes place in the local 
field determined by the second group of dimers, either 
unpolarized J dimers for the BEC transition in Jd2, or 
polarized Jd2 dimers for the BEC transition in J. Un- 
fortunately, the critical fields of clinoclase (Fig. |9| are 
too high to observe such effects using present-day high- 
field facilities. Nevertheless, the search for similar sys- 
tems with structurally different spin dimers should be 
an interesting avenue to explore the high-field physics of 
quantum magnets. 



V. SUMMARY 

In summary, we performed a joint theoretical and ex- 
perimental study of the magnetic behavior of the mineral 
clinoclase. Using density functional theory calculations, 
we evaluate the microscopic model for this compound and 



identify two types of spin dimers with the couplings of 
J c:i 700 K and Jd2 — 300 K. Intuitively, one is tempted 
to ascribe them to two types of structural CU2O6 dimers 
in clinoclase. In fact, only Jd2 pertains to the structural 
dimer D2, while the strong coupling J occurs between 
two corner-sharing Dl dimers. Additional couplings be- 
tween the magnetic dimers reach 150 K, but play a mi- 
nor role in the magnetic behavior. Simulations for the 
DFT-based microscopic magnetic model yield excellent 
agreement with the experimental data. 

The magnetic couplings in clinoclase are not solely 
determined by the Cu-O-Cu angles. The ASO4 side 
groups and the hydrogen atoms also play an important 
role by enhancing or suppressing antiferromagnetic con- 
tributions to short-range couplings Jdi, Jb2^ and Jidi- 
Since no hydrogen positions were available from the ex- 
periment, we determined them by optimizing the crys- 
tal structure within DFT. We have demonstrated that 
the magnetic coupling within Dl is strongly affected by 
the hydrogen atom attached to one of the bridging oxy- 
gens. The out-of-plane H position is responsible for the 
almost canceled exchange coupling Jdi- Our findings put 
forward details of the crystal structure, including incon- 
spicuous and typically overlooked effects like hydrogen 
positions, as an important and even decisive factor in 
the magnetic superexchange. 
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FIG. 10. Rietveld refinement of the room-temperature high-resolution XRD data. Ticks 
show the reflection positions of clinoclase. 



TABLE III. The lattice parameters, atomic positions and atomic displacements parameters Uiso 
(in 10 ~^ A^) refined from the high-resolution powder XRD data on the natural sample of clinoclase. 
Positions of hydrogen atoms are determined by a DFT-based structure optimization within GGA. 
All atoms are in the 4e Wyckoff position of the P2i/c space group. 



: 7.2658 A 



b = 6.4588 A 



12.3929 A 



P = 99.487° 



Atom 

Cul 

Cu2 

Cu3 

As 

01 

02 

03 

04 

05 

06 

07 

Hl^ 

H2^ 

H3^ 



x/a 


y/b 


z/c 


t/iso 


0.7864(2) 


0.1400(2) 


0.3296(1) 


0.74(3) 


0.8141(2) 


0.3802(2) 


0.1272(1) 


0.77(3) 


0.3856(2) 


0.3527(2) 


0.4120(1) 


0.99(4) 


0.3083(1) 


0.1509(2) 


0.1794(1) 


0.85(3) 


0.4090(7) 


0.0712(8) 


0.0714(4) 


0.14(6)" 


0.8387(7) 


0.8447(8) 


0.3681(4) 


0.14(6)" 


0.1802(7) 


0.9408(8) 


0.2108(4) 


0.14(6)" 


0.4705(8) 


0.2147(7) 


0.2813(4) 


0.14(6)" 


0.7777(7) 


0.2039(7) 


0.4780(4) 


0.14(6)" 


0.8047(7) 


0.0957(8) 


0.1747(4) 


0.14(6)" 


0.1871(7) 


0.1705(8) 


0.4157(4) 


0.14(6)" 


0.7493 


0.3515 


0.4807 


- 


0.9346 


0.4680 


0.6787 


- 


0.1511 


0.1648 


0.4864 


- 



The atomic displacements parameters of all oxygen atoms were refined as a single parameter 
Hydrogen positions are obtained from the structure relaxation within GGA 
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TABLE IV. The toi and t hopping parameters (in meV) as function of the out-of-plane angle of H, r. Ue-R- 
AFM exchange contributions J^f^ (in K). The ferromagnetic contribution is defined as J™ = Jdi — J^i^- 
calculated with the CU2O6H5 cluster model and a PBEO hybrid functional. 
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FIG. 11. (Color online) Thermogravimetric analysis data taken in a temperature range of 20-500 °C. The mass loss amounts 
to 7.13%. 




FIG. 12. (Color online) The hydrogen bonds (blue) within the 6c-plane of the clinoclase structure. The Dl dimers are shown 
in orange and the As04 tetrahedra in pink color. 



